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General equations are derived for the spectral density and autocorrelation 
functions of a wave train consisting of sine-wave segments with constant 
amplitude. The frequency of a segment may be either /i or / 2 . At regularly 
spaced intervals the frequency is switched or not switched according to a ran- 
dom choice. This type of wave occurs when a random series of marks and 
spaces is sent by frequency-shift keying. The results fall into two main 
classes — namely, that of discontinuous phase at the transitions, which is 
the typical situation in switching between two independent oscillators; and 
that of continuous phase at the transitions, which is more usually applicable 
when the frequency of a single oscillator is changed. Individual treatment is 
given of the various special cases which arise when integral relationships 
between the marking, spacing, signaling, and shift frequencies exist. No 
restriction is made on the relative magnitudes of the different frequencies 
involved. 

I. INTRODUCTION 

The spectral density function, or power spectrum, of a random se- 
quence of signals defines the distribution of average signal power versus 
frequency. This information is useful in system design because it indi- 
cates the frequency band of most importance, the amount of average 
total power in any frequency interval, and the interference which may 
result in other systems. It does not tell us how much distortion the sig- 
nals suffer when the channel does not pass all the frequencies repre- 
sented, nor does it tell us about important spectral components which 
may be associated with unlikely but possible specific signal sequences. 
Keeping these limitations in mind, we still find the spectral density to 
have merit as a descriptive parameter of the system. 

Another important function is the autocorrelation, which is the time 
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domain analog of the spectral density. Because of the Fourier transform 
relationship between the two functions, either can be used as an auxil- 
iary step in computing the other. The autocorrelation function is useful 
in its own right in signal analysis and can be made the basis of control 
operations. 

The present paper presents results on the spectral density and auto- 
correlation functions associated with binary frequency modulation 
systems. There are two general types of operation, which in terms of 
apparatus may be classified as (a) switching between two oscillators, 
and (b) changing the frequency of a single oscillator. The mathematical 
distinction between these two cases will be considered here as shifting 
the frequency with discontinuous or continuous phase respectively. 

The case of discontinuous phase, which is appropriate when we switch 
between independent marking and spacing oscillators, is the simpler one 
to analyze. We assume that the oscillators deliver equal amplitude and 
that each oscillator preserves its own coherence in time, i.e., that the 
two frequencies are constant. The waveforms in intervals containing 
the marking frequency, say, are segments of a sine wave having the 
marking frequency and extending throughout all time. One would 
intuitively expect, therefore, to find discrete spectral lines at the mark- 
ing and spacing frequencies. Analysis verifies that if mark and space 
signals have independent equal probabilities, each of the two discrete 
components has half the amplitude of the complete FSK wave. 

In addition there is a continuous spectrum consisting of switching 
function spectra centered at the marking and spacing frequencies. Since 
the signal wave is discontinuous at the switching instants, the associated 
voltage spectra fall off only as 1// at high frequencies. This means that 
the spectral density function ultimately falls off as the inverse square of 
the frequency. Degenerate cases arise when there are commensurable 
relationships among the marking, spacing, and signaling frequencies. 
If these special relations are such as to produce continuous phase at the 
transitions, the resulting continuity in the signal wave leads to an ulti- 
mate inverse fourth-power variation of spectral density with frequency. 
If in addition the derivative of the phase is continuous, an inverse sixth 
power is obtained at remote frequencies. 

In the case in which the frequency shifting is done with continuous 
phase, the signal wave is continuous at all times. The spectral density 
function must, therefore, fall off at least as fast as the inverse fourth 
power at frequencies remote from the center of the signal band. This is 
in accordance with the well-known fact that frequency-shift keying with 
continuous phase does not produce as much interference outside the 
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signal band as the discontinuous case. The analysis of the continuous- 
phase case is considerably more difficult, even though the final results 
are of fairly simple form. 

It was found necessary to distinguish between four possible cases. In 
the most general of these, in which there are no degenerate relationships 
among the three frequencies involved, the line spectra completely disap- 
pear and the spectral density function is continuous at all frequencies. 
When the difference between the marking and spacing frequencies is a 
multiple of the signaling frequency, defined as the sum of the number of 
marks and number of spaces per second, line spectral terms appear at 
the marking and spacing frequencies, and the continuous part of the 
spectral density function changes its form. The continuous part of the 
spectrum in this case is found to depend also on whether or not the sum 
of marking and spacing frequencies is a multiple of the signaling rate. 
A curious behavior occurs when the frequency shift is an odd multiple 
of half the signaling rate. Here no discrete components appear in the 
spectrum, but the continuous spectral distribution undergoes a sudden 
change relative to that at infinitesimally close values of frequency shift 
not possessing the critical property. 

Table II given in Section VI lists the various cases together with the 
corresponding equation numbers for the associated spectral densities 
and autocorrelation functions. Section V gives illustrative curves of 
these functions for various relations among the marking, spacing, and 
signaling frequencies. 

II. DISCONTINUOUS PHASE 

Before considering the frequency-shift keying problem of this section, 
it is helpful to develop some general results applicable to random 
switching between two waves. Let y(t) be a given function of time 
which is bounded for all t ^ and such that the limits later encountered 
exist. Let x(t) be a random telegraph wave defined by 

x(t) = x n , nT £t < (n+ 1)7' (1) 

where n. = 0, 1, 2, • • • , and the x a 'a are independent random variables 
which assume the values ±1 with equal probability. We calculate the 
spectral density w v (f) and the autocorrelation ff,.(r) of 

v(t) = x(t)y(t). (2) 

We note that v(t) is generated from the source of y(t) by inserting a 
reversing switch for which a random choice between positions is made 
at instants T apart. 
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The spectral density of v(t) is given by the limit of 2<| S„(f) \ 2 )/NT 
as N —> oo. Here < ) denotes "ensemble average" and 
i»jyr 



S N (f) = [ e~ M V(t) dt co = 2rf 

-E^- ianr f dt <T iut y(nT + 0- 

n=n •'n 



(3) 



71 = •'0 

Since (x n x m ) is 1 if n = m and if n t* m 

2 N ~ l 

JV->oo N 1 71= o | Jq 



w v (f) = lim — £ / dt e~ iwt y{nT + l) 



(4) 



The autocorrelation function R v (t) may be calculated either by 
averaging (v(t)v(t + t)> over all t ^ with t held fixed, or by taking 
the Fourier transform of (4) ; thus 

Rv(t) = I w v (f) cos cor df. (5) 

Both methods show that R v (t) = for | r | ^ 2 1 and 

B„(r) = 7* / d/ lim ^ 2 2/(»^ + V(»r + * + r) (6) 

i •'O AT-.00 iV ;i=0 

for ^ t < T. 

Return now to the frequency-shift keying problem and consider the 
signal wave 

mi(0 nT ^ t < (n + l)T 

u(t) = or 

u,(t) n = 0,1,2, •■■ 

(7) 
u k (t) = A cos ( Uk t + 6 k ) k = 1, 2 

where the choice is made independently and with equal probability for 
each interval of length T. The signaling frequency is co s = 2v/T rad/sec. 
The wave u(t) may be written as 

u(t) = u + (t) + x(t)u-(t) (8) 

where x(l) is defined by (1) and 

u + (t) = i Wl (0 + i« 2 (<) 

(9) 

u_(0 = K(0 - Jm(0. 

Thus, u(t) is the sum of two steady-state cosine waves, given by u+(t), 
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and a random component given by 

v (t) = u(t) - u+(t) = x(t)uM)- (10) 

The random component assumes the values ±w_(0- 

When w_(f) is identified with y(t), the spectral density w v (f) of v(t) 
can be obtained from (4). Since the algebra is rather tedious the pro- 
cedure will be merely sketched. The integral in (4) is now 



fdte-^u-inT + t) =±e- iuT,a 



(11) 
{g(u - wi) exp [jOi + juiT(n + §)]+•••) 

where the braces contain four terms similar to the first and 

, sin (aT/2) , . 

(Aa)= (a/2) ■ (12) 

The cosines in u-(nT -f- t) are expressed in exponential form before 
integrating. 

Multiplying (11) by its conjugate complex gives 16 terms. Substitu- 
tion in (4) gives rise to expressions of the form 

. . ,. . N II, if X = 2mt 

7 ( X ) = lmi 7 ,(X) = | Qj x rea , and ^ 2ym 

where m is an integer. A typical value of X is 2oiT. At this stage w v (f) 
is given by 

8rM, ; (/) = ff 2 (w - co,) + f/ 2 (co + wi) + f/ 2 (co - co 2 ) + </ 2 (co + **) 

+ 2 cos (20i + wi7')y(w - wi)jf(« + oji)7(2ojiT) 
+ 2 cos (202 + o*T)g(u — ui)g(u + on)y(2ui 2 T) 

- 2 cos | 2 + 0i + f - ^ - J T fo(« ~ wi)ff(« + <*) (14) 
+ 9-(w + wi)ff(co - oj 2 )]t(w2T + wiT) — 2 cos 2 - 0i 

+ fe-=-^) r]b(« + <ft)0(* + o> 2 ) 

+ {/(w — ui)g(u> — ut)]y(.tatT — o>iT). 
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Some further reduction leads to the final result 



2TA- 2 w v (f) = 



G(w — o)\) . G(co -f" coi) G(u) — a: 2 ) . G(o) 4~ a.'s) 



+ 



+ 



+ 



(w — coi) 2 (w -j- a>i) 2 (co — co 2 ) 2 (w -f- W2) 2 



2 cos 2ftG(u - wi) ( „,. 2 cos 2ftG(a> - a- 2 ) ,„ „v 
— ^— —5 7(2coii ) -i 2 _ it2 7C2a.> 2 i J 



wr 



a>2 



+ 



- 2 cos «?, + ft) r . gtog 

G(( ° ~ a ' 2) l 7 (a- 2 T + W1 T 

[03 + coi) (a> — W 2 ) J 



W2) 

) - 2 cos (ft - ft) 



(15) 



L(w — wi)(co — 



+ 



G(w + a- 2 ) 



where 



co 2 ) (co + w t ) (w + a> 2 ) 



(7(a) = sin 2 (oT/2) 



T(«!^ — 



Wl 



r) 



(16) 



and 7(2a. 1 7') ) 7(2u> 2 T), 7(^2?' + «iT), y(a>tT - u>iT) are zero except 
when 2oji/w 8 , 2oj 2 /o) s , (a> 2 + «i)/«« , (w2 — wi)/w« , respectively, are 
integers (in which case the corresponding 7's are unity). 

In the special case in which 2lo]/lo s = m and 2u 2 /w 3 = I, with / and 
m integers and I + m even, all of the 7's in (15) are equal to unity and 
the G's are equal to sin" (a>7'/2) if m is even and to cos" (w7'/2) if m is 
odd. In particular, when ft = ft = the following simple result is 
obtained 



w v (f) = 



2A (0:2 — o»i )' 



T(o> 2 - Wl 2 ) 2 (co 2 - a>2 2 ) 2 



X 



. 2 wT 7 
sin — , m even 



2 &T u 

cos — , m odd 



:n) 



This spectral density function varies as 1/to for large co, showing that 
the waveform of the signal is continuous and has a continuous derivative. 
The spectral density function of u(t) as denned by (7) differs from 
that of v(t) by the presence of sinusoidal components of amplitude A/2 
and frequencies ui and a> 2 . That is 

wu(f) = w v (f) + i- «(/ - /1) + 4" *(/ - /*)■ (18) 



It can be shown from (6) and (18) that the autocorrelation functions 
R„(t) and R v (t) of u(t) and v(t) respectively are given by the following 
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set of equations 

A 2 A 2 

R u (t) = Rv(t) + -£- COS WIT + g" COS C02T 

R v (r) =0, I r | > T 
8TA~ 2 R v (t) = (T - r)(cos o>,r + cos W27-) 

Sill 0J1T on /rt m\ SH1 W 2T 



0,08 2^7(20,^) - ^^21 C os20 2 7(2o> 2 T) 



0>l 0>2 



(19) 



, 2(sin W2T + sinojir) /« i a \ / m i t\ 
-}- — =— -^ cos (0 2 + 0i) y{co2i + o)ii ) 

0)2 + 0,1 

. 2 (sin o>2T — sin out) /. a \ ( m t\ f\ < ~ < T 

-| — - cos (02 — d\) 7(oj2i — o»ii ; u s= t =5 i 

0>2 — OJl 

fi„(-r) = R v (t). 
III. CONTINUOUS PHASE 

Consider the signal wave 

A cos (wxt + 0„) »T ^ * < (n + 1)T 
w(0 = or (20) 

A cos (ohjJ + <M n — 0, 1, 2, • • ■ 

where the choice is made independently and with equal probability for 
each interval of length T. The initial values at t = of the phase are 
() = O = <£, and the succeeding values n , <j>„ are to be chosen so as to 
make the phase of u(t) continuous at the transition points. 

Let 

a = i(o, 2 + Wl ) /3 = £(tt*-wi). (21) 

Then 

wi = a - o, 2 = « + 0. (22) 

Set 

u(0 = .1 cos B n (t), nT ^ t < (n + 1)T, n = 0, 1, 2, • • • (23) 

B (t) = (a + .r,/3)£ + <£ (24) 

B„(tf = at + x, l+] 0(t - nT) + + pT £ av , n > 0. (25) 
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We assume xi , x% , • • • to be independent random variables, each of 
which is equally likely to have the value +1 or — 1. We verify that 
within the interval beginning at t = nT, the frequency is a + x n +i&, 
which is equal to 02 if %n+i = +1 and equal to ui if x„+i = — 1. There- 
fore, the function B n {t) satisfies the condition of an equiprobable 
choice between the two frequencies in each interval. We also note that 
the phase at the beginning of the typical interval is 

B n (nT) = anT + + 0T £ x r (26) 

while the phase at the end of the previous interval is 

7.-1 
B n -i(nT) = anT + <f> + x n p[nT - (ft - \)T] + pT £ x T 



(27) 



= B n (nT). 



Thus the function B„(t) also satisfies the required condition of continuous 
phase. 

We have evaluated the spectral density function of u{t) by two 
different methods, namely 

(a) by taking the Fourier transform of the autocorrelation function, 
and 

(b) by direct evaluation from the Fourier transform of the signal 
wave over a long time interval. The two methods are of comparable 
difficulty. The same results are finally obtained by both procedures, 
although the agreement is not immediately evident from the expressions 
which emerge naturally from the two sets of calculations. 

To evaluate the autocorrelation function, we first calculate the 
average value of u(t)u(t + t) over the ensemble at fixed t. Set t ^ 
and define k as the member of the set 0, 1, 2, • • • satisfying the inequality 

(n + k)T ^ t + t < (n + k + l)T (28) 

with n defined by nT ^ t < (ft + \)T. Let E k (t, t) represent the 
mathematical expectation of u(t)u(t -f- t) with / and t fixed. If k = 0, 
the values of t and t + t lie in the same signaling interval and the same 
function B,,(t) applies to both. When k > we use the function 
B n (t) for u(t) and the function B n+k (t) for u(t + t). 
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We calculate 

= 2<cosB w (0 cos B n (t + r)) 

= Re <[exp (jB„(t + r) - jB H (t)) + exp (jB n (t + r) + jB n (t)))) 

= Re (\ exp (jar) exp (> b+ ,0t) + exp (ja(2/. + r) + j2tf>) ^ 

• exp (jx n+i p(2t + r - 2nT)) ft exp (j2.r r /3T)j) . 

Since the x's are independent, the average of a product of functions in 
which the variables appear separately is equal to the product of the 
averages of the individual functions. Since each x has only two possible 
values with a probability of one-half for each, we evaluate the expecta- 
tions of individual terms by inserting the sum of the two possible 
functions with weighting factor one-half. Performing the necessary 
operations, we find 

2 
"Tj E (t, t) = cos «t cos /3t ^0) 

+ cos (2ar -f ar + 20) cos P(2t + r - 2nT) cos"2/37\ 
The corresponding calculation for Ic > leads to the result 

ja E k (t, r) = cos 0{t + t - 7iT - kT) 

cos*- 1 /3!T[C0S ar cos 0(t - nT - T) + COS (2at + ar + 2<(>) (31) 

cos pit - nT + T) cos" 2/3T]. 

The lag interval t is bounded between adjacent multiples of T by 
denning m as the member from the set 0, 1, 2, • • • which satisfies mT ^ 
T < ( w + 1 ) T. We observe that the number k denned by (28) is related 
to m as follows 

nT ^ t < (m + n + \)T - t 
m + 1, (m + n + l)T - r ^ t < (n + l)T. 

The autocorrelation function R„(t) is the average of (u(t)u(i + t)) 
taken from t = to t = «>, i.e. 
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R,M = lim — , / (u(t)u(t + r)> dt 

jv-oo Jy 1 J o 

= Km =£= £ / E m (t,r) dt 

AT-»oo IV I n=0 J nT 



4- / J? m+ i(/,r) A (33) 

J (m+n+l)T-T J 



(ni+l)T-r 



1 n— i i nnti)'-T 

= lim — , £ / #„(/ 4- nT,r) rf/ 



+ / # m+1 (/ + n!T,r) #1 . 

J (m+1)T-T J 



As indicated in (30) and (31), the case of m = 0, i.e., < r < T, 
requires a separate treatment from that of m > 0. In order to perform 
the integrations indicated in (33) it is convenient to expand (30) and 
(31) into the sum of terms in which t appears only once. The expanded 
equations are, with ^ t < T 

2 
■j2 E (t + nT, t) = cos <xt cos fir 

+ § cos" 2/3T cos [(a + (8)(2< + r) + 20 + 2na!T] 

+ i cos" 2(3T cos [(a - j8)(2i + r) + 20 + 2naT] (34) 

4 

-2 E k (t 4- nT 7 , r) = cos ar cos fc_1 (ST [cos $(t + T - fcT) 



(35) 



4- cos 0(2* 4- r - r - fcT)] 

4- cos fc_1 pT cos" 2/3T {cos 0[t - (k + 1)T\ cos [a(2* + t) 

4- 20 4- 2naT] 4" § cos [(a + 0)(2J + r) + 20 4- 2naT - (k - 1)/3T] 

4- \ cos [(a - /3)(2« + t) + 20 4- 2ncxT 4- (A - 1)/MT). 

We note the possibility that some of the terms will not contribute 
anything to the result after the limiting process in (33) is performed. 
The expressions in (34) and (35) can be divided into two classes: those 
which do not contain n and those which contain n in the form of a 
factor of type cos" 2/3!F cos (\p 4" 2naT). In the former group, we sum 
N equal terms and divide by NT; hence the summing and limiting 
operations are equivalent to a division by T. In the other group, the 
limit is zero if the sum remains finite as N becomes indefinitely large. 
The only case in which the sum does not remain finite is that in which 
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the terms are equal for all n. Equality occurs if aT and {ST are both 
multiples of t and also if aT and 0T are both odd multiples of v/2. If 
neither of these conditions exists, the contributions of the terms which 
depend on n are zero. We shall call this the incommensurable case and 
shall treat it first. Afterward, we shall consider the effect of commen- 
surable relations. 

Except for the cases we have specifically excluded, we can now write 
for < r < T 

1 2 pT—r 



A f 

Ru(t) = ^ / COS OCT COS 0T (ti 



(37) 



2T J 
+ ^ ( C os ar[cos 0t 4- cos 0(2/ 4- t - 2T)] dt (36) 

4 2 

= -—=, [0(2T - r) cos 0r + sin 0r] cos ar. 
40 1 

¥ov t > T 

a 2 /-(m+DT-r 

Ru ( T ) = £L / C os ar cos"- 1 0T[cos ${ T + T - mT) 

4i •'o 

+ cos 0(2/ 4- t - T - mT)] d/ 

4 2 C T 
_|_ ^ / cos ar cos"* 0T [cos 0(r - mT) 

4T *'(m+l)T-r 

+ cos 0(2/ 4- r - 2T - mT 1 )] dZ 

= -^ cos ar cos'" -1 (3T[pT COS 0(r 4- T - mT) 

4- 0(r - mT) sin 0T sin 0(r - mT) 
4- sin 0Tcos0(r - mT)]. 
The spectral density function w„(f) is given by 

Wu (f) = 4 f « U ( T ) cos a>r dr. (38) 

It is convenient to evaluate the integral in two parts 

4 f R u (t) COS urdr = ^-= I [0(2T - r) COS 0r 4- sin 0r] 

J o 43 i J o (39) 

■[cos (co 4- a)r 4" COS (a) — a)r] dr 

and 



(40) 
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4 / R u (t) COS cot dr 
J T 

= Ss E / cos «t cos'" -1 pT[pT cos /3(r 4- T - mT) 

pi m=l ■'mT 

4- (8(t - mT) sin pT sin B(t - mT) 
+ sin BT cos /3(t — mT)] cos wt dr 

4' J r 7 ' 

= '— / [pT cos B(t + T) 4- pV sin BT sin pV 
pi •'u 

4- sin pT cos Br] JT cos a(r + mT) cos co(r + //iT) cos'" -1 BT 

m=l 

= 4- f fo^ cos B(t+ T) + 8t sin pT sin p> 

+ sill pT cos Bt][G(u + a, r) + C7(a> - a, t)] dr 

where 

r?f« ) cos y( T + T) - cos pT cos yr , , 

KJ,TJ " 1 4- cos 2 BT - 2 cos BT cos t/7 7 ' ^ ; 

The summation is performed by writing 

2 cos <x(t + mT) cos o>(t 4- mlT) cos"' -1 /3jT 

as the real part of 

[exp (j(« + «)(t + mT)) + exp (j(co - a)(r + mT))] cos"- 1 /3T 

and thereby obtaining a geometric series. The series fails to converge 
when the absolute values of the individual terms are unity. For this 
reason, we must now exclude the case of BT equal to a multiple of ir 
no matter what aT is. We have thus accumulated three special cases of 
commensurability to be given individual treatment later. 

The remainder of the calculation is straightforward, but appears to 
lead into a morass of complication. The key to an end result of pleasing 
simplicity, which the autocorrelation method tends to conceal, is to 
arrange the work as follows 

H (to + ql) 

1 + cos 2 BT - 2 cos BT cos (a + a)T 

, H(a> - a) 1 

1 4- cos 2 BT - 2 cos BT cos (o> - a)T J ' 



«".(/) =f[ 



(43) 
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The integral (40) splits naturally into this form. In the integral (39) 
we associate the term cos (« + a)r with the first part of (42) and the 
term cos (« — a)r with the second part. Then 

H(y) = 1- ( {[14- cos 2 pT - 2 cos pT cos yT] 

2/3 J o 

■[p(2T - T ) cos p T + sin pr] cos yr 

4- [pT cos p( T + D + ^r sin 0T sin /St 

4- sin pT cos /3r] [cos ?/ (r + T) - cos pT cos yr]j <2t 

= 1- ( {(l + cos- 0T)\fi(2T - r) cos /3t 4- sin j8t] 

2/3 ^o 

- cos /3T[/3T cos 0(t + 7') 4- 0t sin 02' sin 0t 
4- sin pT cos /3r] } cos ijt (It 

+ — [ {pT cos /3(t + T) 4- 0t sin /37' sin fir 
2/3 *'o 

+ sin pT cos /3r - cos pT[p(2T - r) cos Pt 
+ sin /8r] J cos y (t + T) dr 

- C -^ff T [P(2T -r) COS fir 

+ sin p T ] cos t/ (t - T) (It. 

In the second integral, substitute r 4- T = / and in the third integral 
substitute T - t = t'. Dropping the primes after the substitution and 
combining terms where possible, we then find that the result can be 
written in the form 

\ C T 1 C' 1T 

H{y) = -±- I hi(r) COS yr (It + —I h,( T ) COS \Jt (It (44) 
2/3 J» J.p j t 

where 

/i,(r) = 2p(T - t) cos Pt - Pt cos /3(t - 27 1 ) 4- sin /3(t - 27') 

4- 2 sin Pt (45) 
A 2 ( r ) = j8( T - 22 1 ) cos/3(t - 2T) - sui/3(t - 22 1 ). (46) 
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The integration in (44) leads to the result 

^..^f^f^-jlJ. (47) 

The complete equation for the spectral density can now be written in 
the form 



w u (f) = 



2A 2 sin 2 h ~ Ml ) T sin 2 ( a ~ " 2 ) T 
T\l - 2 cos (<a - a)T cos /?T + cos 2 /37'| 

1 



■ 



O) — COi CO — 0>2_ 



+ 



2A sni V — 2 — >? 7 ™ V — 2 — >> ? 

7'(1 - 2 cos (w 4- a)T cos /ST 1 + cos 2 07'] 

1 1 



(48) 



CO + CO] CO + C02_ 

The intermediate step converting the original integral (43) to the form 
(44) is very helpful in reducing the labor required to obtain the final 
form (47). Incidentally, (44) shows that the function H(y) is the 
Fourier transform of a function of r which is time-limited to the range 
to 27 1 . We also point out that discrete components do not appear in 
u(t). The spectral density function is continuous at all frequencies and 
varies as the inverse fourth power of the frequency at frequencies remote 
from coi and co 2 . The latter property must exist because the waveform 
of the signal is continuous at all times. 

We now return to the three cases of commensurability which we 
found necessary to avoid in deriving the general result of (48). When 
(IT is a multiple of w, an examination of (25) shows that B„(l) differs 
from at + a*„ + i/3/ + <f> by a multiple of 2ir, and hence 

u(0 = A cos (at + s»+i# -I- <f>), nT ^ t < (n + 1)7'. (49) 

Comparison with (7) shows that we now have one of the degenerate 
cases in which the generally discontinuous phase becomes continuous: 
i.e., that case in which 0i = &•> = <f> and coo — oh = rta, , r being an integer. 
When (3T is a multiple of t and aT is not a multiple of ir, we have 
co2 — cot = rco s , co 2 + coi ^ la) s , and it follows that 2coi/co„ , 2co 2 /co s are 
not integers. Setting t(2co 1 7 7 ), 7(2co 2 T), y(u) 2 T + o^T) to zero in (15) 
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and using 

*- , - - '[( ! T ! ) , ]"* f [(T) ,+ l] 

. 2 \<a — ct\ m 

sin — - — i , r even 

cos 2 r^^l;r j rodd (50) 

C(o} — o> 2 ) = (t(o) — o>i ) (51) 

together with the corresponding expressions for (7(o> + o>i), (?(o> + wo), 
gives «)„(/) in 



Wu(f) = W,(f) + i «(/ " /l) + «(/ " /,)] 



(52) 



where the notation is the same as in (18). When r is an even integer 
we calculate 



, (/) = dl W ( "-<* ) T \-l 1 T 

2T [ \ 2 / |_co — o>i co — C02J 

+ ^- + -Vf » L_T| 

\ 2 / [_u + W] 0) + o> 2 J J 



(53) 



and when r is odd 



27' ( \ 2 / |_C0 — COi 10 — 0)2 J 



+ cos' 



(!L±^ r T_i L_T\ . 



(54) 



In the next case both /3T and aT are multiples of ir; i.e., o> 2 — o>i = rw, , 
C02 + o>i = /o) s , and 2o!i/o s , 2o) 2 /o) s are integers. Now all of the terms in 
(15) must be considered and 



G'(o) — 0)1 ) = G(co — OJ 2 ) 

= (j(o) + 0)1) = r/(o) + W2) = 



. 2 0)T „ . 

sin — tor I — r even 

— 

cos 2 — for Z — r odd. 



(55) 



Combining terms in (15) gives 
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. 2 OiT . , 

sin — - for / - 


■ r even 




4 2 


2 




u;,(/) 


" 2~f 


2 uT c , 
cos — for t - 


- r odd 



[(—--4)' 

(_\C0 — 0>l CO — OJ2/ 



(56) 



+ P-— 4-Y 

\C0 + COl CO + 402/ 

+ 2 C_J Lj) ( ' 1 ) cos 2 J . 

\0> — OJl CO — C0 2 / \C0 + COi CO + C02/ J 

In the last of the exceptional cases, both aT and pT are odd multiples 
of ir/2. That is, 

2aT = (2/+ 1)tt, Z = 0, 1,2, ••• 

2/37' = (2r + 1)tt, r = 0, 1,2, •••. 
Equivalently 

co 2 /co, = (l + r+ l)/2 
coi/co„ = (I - r)/2. 
The value of 7?., obtained by substituting (57) in (.34) is 



(57) 



(58) 



7sBo(« + »7\ t) = cos ar cos 0t + \ cos [(« + 0)(2* + T ) + 20] 
^ (59) 

+ ^cos[(« - (3)(2«+ t) + 20]. 

Substituting /«• = 1 in (35) and then inserting the special conditions of 

(57) we obtain 

4 

72 Exit + nT, t) = cos ar[cos @t - cos 0(2* + t)] 

- cos jSr cos [a(2t + r) + 20] + § cos [(a + 0)(2* + r) + 20] (60) 

+ £cos[(a - j8)(2/+ t) +20]. 

Since I3T is an odd multiple of v/2, the value of cos @T is zero. For 
k > 1 the right-hand member of (35) contains cos /37 1 as a factor. Hence 
E k (t + nT, t) vanishes for k > 1. It follows that the autocorrelation 
function vanishes for r > 2T and there can be no discrete sinusoidal 
components in the spectral density function. 

A better understanding of this remarkable behavior can be obtained 
by examination of a particular signaling interval in which we are equally 
likely to find one of the two possible waves A cos (coif + fa) and 
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A cos (io 2 t + ft). At the next switching instant, say t = nT, we either 
continue with the same wave or shift to the other frequency with 
continuous phase. There are thus four possible waves in the succeeding 
interval, viz., 

(»') A cos (ant + ft) 

(ii) .4 cos [coo/ + ft — (coo — <a\)nT\ 

(Hi) A cos [uit + ft + (coo — coi)nT] 

(iv) A cos (u 2 t -f ft). 

Since (co 2 — ui)T is an odd multiple of ir, the second and third terms 
can be written as (-)".4 cos (coot -f- ft) and (-)"A cos (co^ + ft) 
respectively. Waves (i) and (ii) are possible when the initial frequency 
is coi , and waves (Hi) and (iv) are possible when the initial frequency 
is coo . Now examining the possibilities after the next subsequent switching 
instant, t = (n + 1)T, we find that for the original frequency equal to 
coi , the waves (i) and (ii) can change to the four possible waves: 



(i) A cos (uit + ft) 

(n) (-)" +1 /l COS (ttrf + ft] 

(m) ( — )".4 cos (co 2 t + ft) 

(«•) — A cos (coi/ + ft). 



It will be noted that the first and fourth waves are the same except for 
opposite signs and likewise for the second and third. In other words, for 
any observed value of frequency in a specified interval, the possible 
waveforms in the second succeeding interval can be divided into equally 
likely positive and negative matching pairs. This behavior, once estab- 
lished, must continue into all succeeding intervals. Hence the average 
lag product over the ensemble at fixed t and t must vanish for r greater 
than 2T. 

This may also be seen by noting that when 2f3T = (2r + 1)tt, the 
quantity 0T(xi + • • • + .r„ - nx n+1 ) appearing in the definition (25) 
of B n (t) is an even or odd multiple of ir according to whether 
(•t'i + • • • + x n — nx n+1 )/2 = r„ is even or odd. Hence 

A(-Y n COS (w 2 < + <t>), Zn+l = 1 

«W = ., r (61) 

.4(-) r "cos (co,« +</>), .r„ +1 = -1. 

For any observed set of x x , •■■ x n+1 leading to a definite u(t) in the nth 
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interval, the waveforms in the second succeeding interval can be written 

as 

A(-) r » +(r "+ 2 - r " ) COS M + <*>), Zn+3 = 1 

ii(-) r » +(r "+ 2 - rn) COS («rf + 0), X n+3 = -1. 

Since r„ +2 - r n contains x n+2 only through the term x n+2 /2, the forms 
in (62) are of random sign independent of the original form (61). 
Hence the waveform in the second succeeding interval is entirely inde- 
pendent (in sign and frequency) of the waveform in the original interval. 
Since E (t + nT, t) and Ei(t + nT, t) in (59) and (60) are found 
not to depend on n, we can evaluate the autocorrelation function by 
averaging over t in a single signaling interval as follows 

R»(,t) =i f ^ E (t + 7lT,r)dt, 
1 Jo 

+ A /* E^t + nT, T )dt < r < T (63) 

T J T-r 

Ru { T ) = i f 2 ' " Jfc(l + nT, r) dt, T < r < 2T. (64) 

For < r < T, we calculate 

^ B«(r) = [(2T - r) cos /3r + ^-^ cos ar 

A L * J (65) 

fcos /3t sin ar s in (a + 0)t _ sin (a - /3)t ~| 
+ cos2<^ 2 (a + 0) 2(«-0) J' 

For T < r < 27\ 

«B.(r)-r(2r-r)«.fr + ^l««- 

A L P J ( 66) 

fcos /3t sin ar sin (a + 0)t _ sin (a - fl)r ~l 
_ CO s2<^ - 2 (a + P) 2(a-0) J" 

The spectral density function is then found to be 

/.2T 

WuU) = 4 / Ru(t) COS cot^t 
Jo 

= 2 sinl^Tr/_l L.Y + (-1 L_Y (67) 

2T [_\0) — COl W — 0)2/ \W + CO! 10 + C0 2 / 

+ 2 (_! L_) ( 1 1 ) COS 2 /| . 

\C0 — OJi CO — C0 2 / \« + «1 W -f- ^2/ J 
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The difference between this result and the limit of the general expres- 
sion (48) when the special case (57) is substituted consists of the term 
containing 2<£. This exceptional case thus has the property of remem- 
bering the initial phase angle even though the spectral density is con- 
tinuous. The reason is that the phase of the wave with respect to the 
signaling interval must remain fixed throughout all time and is, there- 
fore, not subject to averaging as in the more general case. 

IV. FOURIER TRANSFORM METHOD 

Expressions (48) and (67) for w u (f) have been obtained by first 
computing R u (t) and then taking its Fourier transform. As mentioned 
earlier, w u (f) can also be obtained by working with a Fourier-type 
integral of u{t) taken over a long time interval. This method will now 
be sketched for the case in which 2/3 T is not a multiple of jr. Most of 
the intermediate steps are omitted. If they were included, the length 
of the derivation would be comparable to the derivation based on 
Ru(r). 

The spectral density w u (f) of u(t) = A cos B n (t) is the limit of 
2<| S N (f, NT) \ 2 )/NT as N -> « . In this expression 



(08) 



S„(f,NT) = f c- Jut u(t)dt=^2 Sn 

s„ = Ae~ ianT I e~ jut cos B n (t + nT) dt (69) 

•'0 

where s„ is a function of /. The ensemble average of 
I S K (f, NT) | 2 = Z Z s n s m * 

n=0 m=0 
N-l JV-1 N-k-1 

= E S»8»* + Z Z (S n s n+k * + S n+k S n *) 
n=0 A=l n=0 

is the sum of terms of the form 



(70) 



(s n+k s,*) = A 2 e' iukT f T e- jut dt [' e juT 

(cos B n+k (t + nT + ¥T) cos B n ( T + nT)) dr. 



(71) 



In these equations s,* denotes the conjugate complex of s n . 

The procedure used to obtain expressions (30) and (31) for E (t, t) 
and E k (t, t) leads to 
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2(cos B n (t + nT) cos B n {r + nT)) = cos a(t - r) cos 0(1 - r) 

+ cos [a(i + r + 2»T) + 20] cos 0(J + r) cos"2/5T 

2<coa B.^(( + »T + tr) cos B„(t + nT)) = cos «(/ - r + AT) 

cos 0t cos 0(t - r) cos'" 1 0T + cos [a(J + r + 2»T + AT) (73) 

+ 2<t>] cos fit cos /3(r + r) cos* -1 0T cos" 20T 

where & > in the last equation. 

We shall consider only the case in which 2/3T is not a multiple of tr. 
Then the terms in (72) and (73) containing cos" 2/3T contribute nothing 
to the left-hand sides of 



(74) 



(75) 



lim iV _1 £ (s n s,*) 

W-»oo n =0 

= 2~ l A 2 / e~ iut dt \ e 3UT dr cos a(t - r) cos fi(t - t) 

Jo •'0 

N-k-1 
lim iV" 1 E (Sn + kS n *) 

= 2- 1 AV i "* r cos fc ~ 1 ^r [ e~ M dt [ e iuT cos a(t - t + AT) 

cos /3i cos j8(t - T) dr. 

Expression (48) for the spectral density w u (f) is now obtained by 
performing the integrations and then summing with respect to /c as 
indicated in (70). 

V. ILLUSTRATIVE CURVES 

Fig. 1 shows a typical curve for the spectral density function when 
the phase is discontinuous at the instants of transition. The curve is 
calculated from 

(rUA % S(f-fi) | &(f-f*) 
w u (f)/A = 5 H £ 

8 8 n (76) 

+ 2T(<a - «j) a 27 , (o J - o> 2 ) 2 

which is an approximation to (15) and (18). It holds when w 2 - «i is 
not a multiple of u, , and in addition a* and w 2 are so large that the 
portion of the spectrum folded back from u = 0, i.e., the portion de- 
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FREQUENCY RATIO, I-I, = (f-f,)/f s 

Fig. 1 — Spectral density of random binary FSK wave with discontinuous 
phase at transitions. Frequency shift = 0.8 times signaling frequency. 



pending on inverse powers of « + a>i and u> -f- o; 2 , is negligible. The 
contribution of the neglected terms becomes appreciable only when the 
marking or spacing frequency is less than the signaling frequency. It is 
convenient to let x = a>/a>* = ///, , x\ = fi/f, , and x-< = f»/f, . 

Fig. 1 is calculated for the case co 2 — oji = 0.8cj 8 , i.e., x% — X\ = 0.8. 
The abscissa is x — x\ . The ordinate is f 8 w u (f)/A 2 . The curve is sym- 
metrical about x — .ri = 0.4. The steady-state terms are represented 
by spikes of infinite height and infinitesimal width at x = Xi and x = x% . 
Each of these spikes has an area of J. The area under the continuous 
curve is j. The total area is \, which is the mean-square value of the 
signal wave per unit of squared amplitude. 

Fig. 2 shows a case of continuous phase corresponding to a frequency 
shift equal to 0.8 times the signaling frequency. This curve was computed 
from (48) where, again, the terms containing inverse powers of <a + wj 
and w + o> 2 are assumed to be negligibly small. Since the infinite spikes 
representing steady-state components are absent, the total area under 
the curve must be \. We note that peaks of finite height and width 
appear just outside the interval bounded by the marking and spacing 
frequencies. These peaks become more pronounced and move toward 
the marking and spacing frequencies as we approach the commensurable 
case in which the frequency shift/.. — j x is exactly equal to the signaling 
rate f s . Fig. 3 shows the curve for frequency shift equal to 0.95 times 
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Fig. 2 — Spectral density of random binary FSK wave with continuous phase 
at transitions. Frequency shift = 0.8 times signaling frequency. 

the signaling rate. Here the peaks are almost twenty times as high as 
in Fig. 2. The limiting case of .r 2 — x% = 1 is exhibited in Fig. 4. The 
finite spikes of Fig. 2 and 3 have now become full-fledged impulses of 
infinite height, infinitesimal width, and area f. They represent the 
mean-square value of steady-state components at the marking and 
spacing frequencies. The continuous part of the curve was calculated 
from (54), noting that (« - a)T = 2ir(x - x x - \) in this case. Fig. 
5 shows a representative curve on the other side of the limiting case, 
with the frequency shift taken equal to 1.2 times the signaling rate. 
The finite peaks now appear inside the interval between marking and 
spacing frequencies. 

It is instructive to study the transition from Figs. 2 and 3 to Fig. 4. 
Since the curves are symmetrical about x — X\ = (.r 2 — X\)/2, it is 
sufficient to consider the region of rapid change near x = X\ . Setting 
Xj - xi = x d , we approximate (48) for w u (f) in this region by 



Wu(f) 



A 2 sin 2 (x — x\)tt sin 2 (x — x 2 )w 

2ir' 2 f s (x — Bi) 2 [l — 2 COS (2.T- — Xn — Xi)lT COS X d T 

+ COS 2 XdT] 



(77) 
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Fig. 3 — Spectral density of random binary FSK wave with continuous phase 
at transitions. Frequency shift = 0.95 times signaling frequency. 
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Fig. 4 — Spectral density of random binary FSK wave with continuous phase 
at transitions. Frequency shift = signaling frequency. 

2377 



2378 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1963 




-0.8 



-0.4 0.4 0.6 1.2 

FREQUENCY RATIO, I-I, = (f-fi)/f s 



Fig. 5 — Spectral density of random binary FSK wave with continuous phase 
at transitions. Frequency shift =1.2 times signaling frequency. 

We are interested in the behavior of (77) as x 2 - Xi approaches unity. 
Setting x% — Xi - 1 - e and x - xi = y, we find that when both e 
and y are small compared with unity we can approximate (77) by 

(y + 2 



Y = f,w u (f)/A 2 tt- 



(78) 



8 (7re 2 /4) 2 + ill + e/2) 2 ' 

It is seen that Y depends on y approximately as shown in Table I. 
The value y = ± °° corresponds to several positive or negative multiples 
of e, and therefore actually becomes small in absolute value as c ap- 
proaches zero. It follows that Y hitches onto the value 0.125 shown at 
.c = .ri in l (, ig. 4. The curves shown in Figs. 2 and 3 correspond to e = 0.2 
and e = 0.05. They show the behavior indicated by (78). In particular, 

Table I — Approximate Ordinates of Spectral Density 
in the Neighborhood of Peak 
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-e/2 
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l/(2irV) 
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J 
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Y obtains its peak value of approximately l/(2w 2 e 2 ) near y = — e/2 
and drops down to about half the peak value at y = — t/2 ± 7re /4. 
When e = 0.05, the peak value of Y is 20.3. The area under the peak, 
as measured by the integral of Y taken from y = — e to y = -\- e, ap- 
proaches | as e approaches zero. This agrees in the limit with the area 
of the impulses shown in Fig. 4. 

The work of Sunde 1 has indicated that the special case in which the 
frequency shift is equal to the bit rate has a theoretical advantage in 
that intersymbol interference can be suppressed at the sampling instants 
in the output of an ideal frequency detector. The results presented here 
show one method by which such a frequency lock can be attained. Since 
the two principal peaks of the spectral density function reach maximum 
height when the condition x-> — X\ = 1 is attained, the output of a 
spectral analyzer can be used to determine the proper bias on the 
tuning control of the keyed oscillator. Another possible instrumentation 
can be devised by use of the autocorrelation function. When the signal 
wave with continuous phase transitions is multiplied by itself delayed 
by a large multiple of the bit interval, the average value of the product 
tends toward zero except when the frequency shift is locked to the bit 
rate. In practice, the advantage of a rigid lock-in to the theoretical 
optimum has not proved to be very significant. The actual reduction 
of intersymbol interference which could be achieved by choosing the 
best value of frequency shift would typically be masked by other 
departures from the ideal conditions. 

Fig. 6 illustrates the case in which aT and $T are odd multiples of 
7r/2. The significantly different properties exhibited are not very pro- 
nounced except when marking and spacing frequencies are sufficiently 
low to be comparable with the signaling rate. The case shown in Fig. 6 
applies when the marking frequency is half the signaling frequency and 
the spacing frequency is equal to the signaling frequency. The frequency 
shift is half the signaling frequency. The curves are calculated from (67) 
for three different values of the initial phase angle. Since there are no 
steady-state components, the area under each curve must be 0.5. The 
peak of the spectral density function changes from 0.763 to 0.857 as the 
cosine of the initial phase angle is varied from — 1 to + 1 . The ordinates 
become zero at x = f, 2, f, • • • . As x approaches infinity the maximum 
values of the intervening loops decrease as x~ when cos <f> = 1 and as 
X for other values of cos 0. 

When the values f\ = f s /2, / 2 = /, corresponding to Fig. 6 are substi- 
tuted in the general equation (48) for w„(f), the result is the case 
cos 20 = shown in Fig. 6. This is to be expected since when the con- 
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Fig. 6 — Spectral density of random binary FSK wave with continuous phase 
at transitions. Frequency shift = 0.5 times signaling frequency. Marking fre- 
quency = signaling frequency. Spacing frequency = 0.5 times signaling frequency. 

ditions/i = /,/2, h = /, are almost (but not quite) satisfied, the phase 
of u = A cos B n (t) changes by a small amount, or by ic plus a small 
amount, from one transition point to another. Over a long period of 
time these small changes accumulate and have the same effect as 
replacing cos 2<j> by its average value 0. 

Figs. 7-12 inclusive show the normalized autocorrelation functions 
corresponding to the cases of Figs. 1-6 respectively. To avoid making 
the curves depend on the values of marking and spacing frequencies, 
we have indicated the envelope of the high-frequency oscillations in 
Figs. 7-11. The autocorrelations are obtained by multiplying the solid 
line curves by %A 2 /cos [(o> 2 + o>i)t/2], which in terms of the lag time t 
is a cosine wave at the midband frequency. The resulting oscillation has 
the value unity at t = and is contained within the solid and dashed 
curves. Fig. 12, which is drawn for specified marking and spacing fre- 
quencies, shows an actual autocorrelation function. 

The typical case of discontinuous phase, which is illustrated in Fig. 7, 
has a linearly damped envelope until t reaches the value T. At time T 
the envelope changes continuously to that of the sum of two cosine 
waves at the marking and spacing frequencies. The latter envelope is a 
cosine wave at half the difference frequency and it persists with un- 
diminished amplitude throughout all values of r greater than T. .Fig. 8 
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Fig. 7 — Envelope of autocorrelation function when phase is discontinuous at 
transitions and f% — f\ m 0.8/, . Actual autocorrelation is full line curve multi- 
plied by (A 2 /2) cos [(w 2 + m)r/2]. 

represents the same case as Fig. 7 except that the phase is continuous. 
The effect is that the envelope of the autocorrelation in Fig. 8 decays to 
zero at infinite lag time instead of oscillating with constant amplitude. 
The decay in each multiple of T after the second one is produced by a 
multiplication of the corresponding values in the preceding interval by 
cos (#2 — Xi)ir, which has the value —0.809 in Fig. 8. As .r 2 — .Ti ap- 
proaches unity, the multiplying factor produces only a slight reduction 
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Fig. 8 — Envelope of autocorrelation function when phase is continuous at 
transitions and fi — fi = 0.8/, . Actual autocorrelation is full line curve multi- 
plied by (A*/2) cos [(« 2 + »i) T /2]. 
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Fig. 9 — Envelope of autocorrelation function when phase is continuous at 
transitions and / 2 — /i = 0.95/. . Actual autocorrelation is full line curve multi- 
plied by (A 2 /2) cos [(a.2 + wi)t/2]. 

in each interval and the oscillations retain appreciable amplitude for 
very large lag times. Such behavior is emphasized in Fig. 9 for the case 
of xi — xi = 0.95, cos (x 2 - xi)tt = -0.9877. The very slow departure 
from constant amplitude oscillations indicates that the signal wave 
contains components which are very nearly sinusoidal. The time domain 
analysis thus agrees with the sharp high peaks found in the frequency 
domain analysis, as shown in Fig. 3 for x 2 — x% = 0.95. 

Fig. 10 shows the limiting case in which the phase is continuous and 
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Fig. 10 — Envelope of autocorrelation function when phase is continuous at 
transitions and / 2 — /i = /» . Actual autocorrelation is full line curve multiplied 
by (A*/2) cos [(« 2 + m)r/2]. 
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Fig. 11 — Envelope of autocorrelation function when phase is continuous at 
transitions and / 2 - /i = 1-2/. . Actual autocorrelation is full line curve multi- 
plied by (A 2 /2) cos [(«> + «i)t/2]. 

Bj _ xi = 1. The appearance of the line spectral terms is indicated by 
the constancy of the amplitude of oscillations for t > T. Fig. 11 for 
x 2 - xi = 1.2 corresponds to Fig. 5. The decay rate is the same as in 
Fig. 8, since cos 1.2x = cos 0.8tt. The period of the oscillations is de- 
creased. 

Fig. 12 shows the singular case in which the sum and difference fre- 
quencies are both odd multiples of half the signaling frequency. The 
values chosen are the same as those of Fig. 6. The autocorrelation func- 
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Fig. 12 — Normalized autocorrelation function when phase is continuous, 
/i = /./2, /« = /., and /* - /i - /./2. 
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tion is time limited, vanishing for all values of t greater than 27\ The 
dependence on initial phase is indicated by the three curves, which are 
drawn for the cases of cos 2<f> = 0, 1, and —1 as in Fig. 6. The total 
variation in the height of the negative peak with <f> is 0.07. 

Fig. 7 differs from Figs. 8-12 in that the discontinuity in phase of the 
signal wave produces a discontinuity in the slope of the autocorrelation 
curve at t = T. This is consistent with the decay of the spectral density 
of Fig. 1 with the inverse square of frequency at high frequencies. The 
spectral densities of Figs. 2 to 6 vary ultimately as the inverse fourth 
power of frequency, which requires not only the slope but the second 
derivative of the corresponding autocorrelation functions, Figs. 8 to 12, 
to be continuous at all values of t. 

VI. SUMMARY OF RESULTS FOR SPECTRAL DENSITY AND AUTOCORRELA- 
TION 

Table II lists the equation numbers of the expressions giving w u (f) 
and Ru(t) for the various cases which can arise. Let/, = 1/T be the 
signaling frequency and /i , fz be the marking and spacing frequencies. 
Also, let I, r denote integers. 

Table II — List of Equations for Spectral Density and 
Autocorrelation of FSK Wave 







Equation 


Numbers 


Case 










VuU) 


R»(t) 


Discontinuous phase: 








(a) general case 




(15), (18) 


(19) 


(b) degenerate cases 




(17), (18) 


(19) 


Continuous phase: 








(c)/j-/x- rf.,f 2 +U* 


If. 


(52), (53), (54) 


(19) 


(d)/2-/. = rf,,f 2 +f 1 = 


V. 


(52), (56) 


(19) 


(e) /. - /, - (r + i)f.) 




(67) 


(65), (66) 


(f) all other continuous phase cases 


(48) 


(36), (37) 



VII. other related publications 

Jenks and Hannon 2 have given spectral density curves for the case of 
frequency shift equal to bit rate which they state have been taken from 
a forthcoming paper by Pushman in the Journal of the British Institute 
of Radio Engineers. The curves shown are in agreement with ours for 
the same case. We have not seen the complete work. Our interest in 
the problem was initially stimulated by discussions with I. Dorros, 
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who has made use of some of our results in a study 3 of the transmission 
of binary data by FM over a band-limited channel. Since completing 
the work, we have become aware of a publication by Postl, 4 who has 
calculated the spectral density for binary continuous phase narrow-band 
FSK in which the midband frequency is large compared with both the 
frequency shift and the signaling rate. 
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